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Abstract 

Berkeley conducts 7 SET! programs at IR, visible and radio wavelengths. Here we review two of the newest efforts, 
Astropulse and Fly's Eye. 

A variety of possible sources of microsecond to millisecond radio pulses have been suggested in the last sev- 
eral decades, among them such exotic events as evaporating primordial black holes, hyper-flares from neutron stars, 
emissions from cosmic strings or perhaps extraterrestrial civilizations, but to-date few searches have been conducted 
capable of detecting them. The recent announcement by Lorimer et al. of the detection of a powerful (» 30 Jy) and 
highly dispersed (« 375 cm"^ pc) radio pulse in Parkes multi-beam survey data has fueled additional interest in such 
phenomena. 

We are carrying out two searches in hopes of finding and characterizing these jus to ms time scale dispersed radio 
pulses. These two observing programs are orthogonal in search space; the Allen Telescope Array's (ATA) "Fly's 
Eye" experiment observes a 100 square degree field by pointing each 6m ATA antenna in a different direction; by 
contrast, the Astropulse sky survey at Arecibo is extremely sensitive but has 1/3,000 of the instantaneous sky coverage. 
Astropulse's multibeam data is transferred via the internet to the computers of millions of volunteers. These computers 
perform a coherent de-dispersion analysis faster than the fastest available supercomputers and allow us to resolve 
pulses as short as 400 ns. Overall, the Astropulse survey will be 30 times more sensitive than the best previous 
searches. Analysis of results from Astropulse is at a very early stage. 

The Fly's Eye was successfully installed at the ATA in December of 2007, and to-date approximately 450 hours of 
observation has been performed. We have detected three pulsars (B0329H-54, B0355-I-54, B0950-H08) and six giant 
pulses from the Crab pulsar in our diagnostic pointing data. We have not yet detected any other convincing bursts of 
astronomical origin in our survey data. 
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1. Introduction 

The Berkeley SETI group conducts seven searches at 
visible, IR and radio wavelengths, covering a wide vari- 
ety of signal types and spanning a large range of time 
scales: SETI@home searches for radio signals with 
time scales ranging from mS to seconds. SEVENDIP 
searches for nS time scale pulses at visible wavelengths. 
ASTROPULSE and Fly's Eye search for dispersed uS 
and mS time scale radio pulses from extraterrestrial 
civihzations, pulsars, or evaporating primordial black 
holes. SERENDIP and SPOCK search for continuous 
narrow band signals in the radio and optical bands re- 



spectively. DYSON searches for infrared excess from 
advanced civihzations that use a lot of energy. 

SETI@home H, ASTROPULSE and SERENDIP V 
are radio sky surveys at the 300 meter Arecibo tele- 
scope. Commensal observations have been conducted 
almost continuously for the past ten years and are ongo- 
ing. Most beams on the sky visible to the Arecibo tele- 
scope have been observed four or more times. We rank 
SERENDIP and SETI@home candidate signals based 
on the number of independent observations, the strength 
of the signals, the closeness of the signals in frequency 
and sky position, and the proximity to stars, planetary 
systems, galaxies, and other interesting astronomical 
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objects. SETI@home uses the CPU power of volun- 
teered PCs to analyze data. Five million people in 226 
countries have participated. Combined, their PCs form 
Earth's second most powerful supercomputer, averaging 
482 TeraFLOPs and contributing over two million years 
of CPU time. Here we describe two of the newest SETI 
searches, Astropulse and Fly's Eye. 

2. SETI Pulse Searches 

One common assumption of SETI is that an alien 
civilization wishing to make contact with others would 
broadcast a signal that is easily detected and easily dis- 
tinguished from natural sources of radio emission. One 
way of achieving these goals is to send a narrow band 
signal. By concentrating the signal power in a very nar- 
row frequency band, the signal can be made to stand out 
among the natural, broad-band sources of noise. By the 
same token, signals leaked from a civilization's internal 
communications may be narrow band, but would be sig- 
nificantly weaker than a direct attempt at extraterrestrial 
contact. 

Because of this, radio SETI efforts have concentrated 
on detecting narrow band signals. When searching for 
narrow band signals it is best to use a narrow search 
window (or channel) around a given frequency. The 
wider the channel, the more broadband noise is included 
in addition to any signal. This broadband noise limits 
the sensitivity of the system. Early systems used analog 
technology to create narrow bandpass filters that could 
observe at a single frequency channel. More recent sys- 
tems use massive banks of dedicated fast Fourier trans- 
form (EFT) processors to separate incoming signals into 
up to a billion channels, each of width ~1 Hz. 

There are, however, limitations to this technique. One 
limitation is that extraterrestrial signals are unlikely to 
be stable in frequency due to accelerations of the trans- 
mitter and receiver For example, a receiver listening for 
signals at 1 .4 GHz located on the surface of the earth un- 
dergoes acceleration of up to 3.4 cm/s^ due to the earth's 
rotation. This corresponds to a Doppler drift rate of 0. 16 
Hz/s. If uncorrected, an alien transmission would drift 
out of a 1 Hz channel in about 6 seconds, effectively lim- 
iting the maximum integration time to 6 seconds. Be- 
cause of the inverse relationship between maximum fre- 
quency resolution and integration time (Av - j;) there 
is an effective limit to the frequency resolution that can 
be obtained without correcting the received signal for 
this effect. 

There are other parameters of the signal that are un- 
known, for example: At what frequency will it be trans- 
mitted? What is the bandwidth of the signal? Will the 



signal be pulsed, if so at what period? Fully investigat- 
ing a wide range of these parameters requires innova- 
tive instruments and enormous computing power. An- 
other possibility that, until now, has not been included in 
SETI searches is that rather than directing large amounts 
of power into a narrow frequency band, an extraterres- 
trial intelligence might direct large amounts of power 
into a narrow time window by sending a short duration 
wide-band pulse. A wide band signal has the advan- 
tage of reducing the importance of the choice of obser- 
vation frequency, and the unavoidable interstellar dis- 
persion of a broad-band pulse reduces the confusion be- 
tween terrestrial and extraterrestrial pulses. In an ion- 
ized medium, high frequencies propagate slightly faster 
than low frequencies at radio frequencies. Therefore RF 
pulses will be dispersed in frequency during their travel 
through the interstellar medium. The dispersion across 
a bandwidth Av is given by 



Av( MHz ) 
^f^(8.3//s) 3 ' dm 
v^'i GHz) 



(1) 



where dispersion measure DM is defined as rigdl and 
is usually quoted in units of cm"^ pc. 

It is possible that a civilization intentionally creating 
a beacon for extraterrestrial astronomers would choose 
to create "pulses" which have a negative DM. Natu- 
ral dispersion always causes higher frequency compo- 
nents to arrive first. A signal in which the low frequen- 
cies arrive first would stand out as obviously artificial. 
As a check on this, as well as to establish our back- 
ground noise limit, we examine both positive and nega- 
tive dispersion cases. Unfortunately, sensitive detection 
of broad-band pulses at an unknown dispersion measure 
also requires enormous computing power 

We are currently conducting two searches for these 
dispersed radio pulses. Astropulse at Arecibo Observa- 
tory and Fly's Eye at the Allen Telescope Array (ATA). 
The Allen Telescope Array is a joint project of the SETI 
Institute and the University of California, Berkeley |3]. 

2.1. Pulses from Extraterrestrial Intelligence 

Astropulse and Fly's Eye consider the possibility 
that extraterrestrials communicate using high intensity 
and wideband but short timescale pulses. It turns out 
that both the conventional (narrowband) method of data 
transfer assumed by SETI@home and the wideband, 
short timescale method assumed by Astropulse and 
Fly's Eye require the same amount of energy to send 
a message. 

In each case, the energy required to send one bit of 
information is proportional to the energy (per area) re- 
quired to send the minimum detectable signal. And this 
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energy is the same for the two methods, so energy con- 
siderations cannot rule out short timescale pulses as a 
medium for extraterrestrials' communication. 

SETI@home's sensitivity in searching for narrow- 
band signals is given by: 



acrTsys 



(2) 



where o; = 2 is a loss from 1 -bitting the data twice, 
cr = 24 is the threshold above noise, Tsys = 28 K is 
the system temperature, G - IQ K Jy ' is the gain, and 
A^poi = 1 is the number of polarizations. 

SETI@home searches for narrowband signals with a 
bandwidth of 0.075 Hz, and its longest integration time 
is (0.075 Hz)"' = 13.4 s. To compute the energy (per 
area) required to send one bit, multiply by Btiai = 1 to 
get 134- \Q--^im-^. 

On the other hand, Astropulse and Fly's Eye are 
searching for short pulses, so the sensitivity of these ex- 
periments will be given in Jansky microseconds, by a 
slightly different formula: 



(3) 



where the variables have the meanings as before, ex- 
cept that now fjnt is the timescale of the pulse and B the 
entire bandwidth. 

For Astropulse, fint is 0.4 and B - 2.5 MHz. Also, 
a - \A because Astropulse 1-bits the data only once, 
A^poi - 2, and cr - 21.5. The resulting sensitivity is 
24 Jy yus. To compute the energy per area required to 
send one bit, multiply by B = 2.5-10^ Hz to get60-10"^*' 
Jm-2. 

For both narrow frequency and broadband pulse 
cases, the expression for the minimum detectable en- 
ergy has the form 



(4) 



with V^fint = 1, SO similar answers are expected. 

Because these energies are comparable, we have 
a chance to detect ETI communications in this new 
regime. 

2.2. Pulses from Evaporating Black Holes 

Another intriguing possible source of short duration 
pulses comes from a suggestion by Martin Rees in 
1977 ISl that primordial black holes, evaporating via 
the Hawking Process, could emit a large electromag- 
netic signature. According to Hawking Q, a black hole 



of mass M emits radiation like a blackbody with a tem- 
perature Tbh given by: 
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This radiation emanates from the black hole event hori- 
zon and comes completely from the black hole's mass. 
For a non-accreting black hole, Stefan-Boltzmann 
yields a lifetime of: 



Tbh - IQ^^^years 



M 



(6) 



For stellar mass black holes, this theory predicts life- 
times of order 10^"* years, much too long to ever expect 
to observe. However, some cosmologies predict the cre- 
ation of numerous small (M ~ lO'^g) primordial black 
holes in the early universe, which according to theory 
could be evaporating now jSj. 

The specific mechanism by which the evaporating 
black hole produces a strong radio pulse has not been 
fully elucidated, but in short, it is thought that the pro- 
cess is similar to the EMP that accompanies super- 
nova explosions. In such a process, a highly conduc- 
tive plasma fireball expanding into an ambient magnetic 
field can exclude the field and create an electromagnetic 
pulse. For typical values of the interstellar magnetic 
field, this pulse would be peaked near IGHz [l]. 

An observation of these pulses would not only pro- 
vide a significant confirmation of Hawking radiation, 
but would also give strong evidence of the existence of 
primordial black holes. 

3. Fly's Eye: Searching for Bright Pulses with the 
ATA 

The Allen Telescope Array has several advantages 
over other telescopes worldwide for performing tran- 
sient searches, particularly when the search is for bright 
pulses. The ATA has 42 independently-steerable dishes, 
each 6m in diameter. The beam size for individual 
ATA dishes is considerably larger than that for most 
other telescopes, such as VLA, NRAO Green Bank, 
Parkes, Arecibo, Westerbork and Effelsberg. This 
means that the ATA can instantaneously observe a far 
larger portion of the sky than is possible with other 
telescopes. Conversely, when using the ATA dishes 
independently, the sensitivity of the ATA is far lower 
than that of other telescopes. 

The Fly's Eye instrument was purpose built to search 
for bright radio pulses of millisecond duration at the 
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ATA. The instrument consists of 44 independent spec- 
trometers using 1 1 CASPER IBOBs. Each spectrom- 
eter processes a bandwidth of 210MHz, and produces 
a 128-channel power spectrum at a rate of 1600Hz 
(i.e. 1600 spectra are outputted by each spectrometer 
per second). Therefore each spectrum represents time 
domain data of length l/1600Hz=0.000625s=0.625ms, 
and hence pulses as short as 0.625ms can be resolvecQ 
We have to-date performed roughly 400 hours of ob- 
serving with the Fly's Eye. Figures [T] and |2] show the 
beam pattern and sky coverage using all 42 antennas. 
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42 antennas; 197.6 square degrees at half power 
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Figure 2: From (2)- The sky coverage of the ATA for an observing pe- 
riod of 24 hours. Both the coverages for southern and northern point- 
ings (corresponding to the respective contiguous regions) are shown. 



Figure 1: From fT\. The beam pattern of the 42 beams at ATA, with 
the diameters equal to the half-power width. This hexagonal pack- 
ing is pointing north. A south pointing results in poor interference 
properties, due to the highly populated areas south of the ATA. 

3.1. System Architecture 

The overall architecture of the Fly's Eye system is 
shown in Figure [3] Each IBOB can digitize four ana- 
logue signals, and 1 1 IBOBs are provided so that 44 sig- 
nals can be processed. The ATA has 42 antennas, each 
with two polarization outputs. A selectiorj^of 44 of the 
available 84 signals is made, and these are connected 
to the 44 iADC inputs. The IBOBs are connected to a 
control computer and a storage computer via a standard 
Ethernet switch. 




12TB RAID Storage 
Server 



3.2. Fly's Eye Offline Processing 

The analysis required for the Fly's Eye experiment 
is, in principle, fairly simple - we wish to search over a 



'Pulses of duration <0. 625ms can also be detected provided that 
they are sufficiently bright, but their length cannot be determined with 
a precision greater than the single spectrum length. 

-The selection is made with consideration for the goal of maximis- 
ing field-of-view - in practice we selected at least one polarization 
signal from every functioning antenna. 




Figure 3: Fly's Eye System Architecture. A selection of 44 analogue 
signals from 42 dual polarization antennas are connected to 44 inde- 
pendent spectrometers implemented in 11 IBOBs. 
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Figure 4: Fly's Eye Rack at the ATA. Two 6U CompactPCI crates (top 
and bottom) house the 1 1 IBOBs. The switch, data recorder computer 
and storage server are all visible. 



wide range of dispersion measures to find large individ- 
ual pulses. Specifically our processing requires that all 
the data be dedispersed with dispersion measures rang- 
ing from 50 cm"^ pc to 2000 cm"^^ pc. At each disper- 
sion measure the data needs to be searched for 'bright' 
pulses. 

The processing chain is in practice significantly 
more complicated than this description suggests. 
Processing is performed on compute clusters, with 
input data formatted, divided and assigned to worker 
nodes for processing. In the worker node flow, the 
data is equalized, RFI rejection is performed, and 
finally a pulse search is performed through the range 
of dispersion measures. The results are written to a 
database where they can be subsequently queried. The 
key feature of the results is a table that lists, in order of 
decreasing significance, the pulses that were found and 
the dispersion measures they were located at. 

Average power equalization is performed on the fre- 
quency spectrum equalized values P'^it). We compute 
the average power over all frequency channels for a 
single integration (time sample t). The power average 
is defined as W(t) = ^ ^ilo' N is the number 

of channels (for Fly's Eye this is always 128). The 
motivation for why it is possible to normalize the power 
is that we expect pulses to be dispersed over many 
time samples, so this procedure should not remove 
extraterrestrial pulses. 

Our strategy for mitigating constant narrowband RFI 
is simply to identify the channels that are affected, and 
to exclude them from further processing. This channel 
rejection is typically performed manually by looking 
at a set of spectra and identifying obviously infected 
channels, which are then automatically excluded in 
subsequent processing runs. 

Intermittent RFI is often quite difficult to automat- 
ically distinguish from genuine astronomical pulses, 
and we followed a conservative approach to try to 
ensure that we do not accidentally excise dispersed 
pulses. Our statistic for intermittent RFI is the variance 
of a single channel over a 10 minute data chunk, 

cr^ = (i 2:^"o ' {P'/(t)f) - ' {PnO)f- Curve- 

fitting determines a cr? outside which it is likely that 
channel / contains time-varying RFI. Future reprocess- 
ing will likely use a more robust method, such as that 
based on a kurtosis estimator ITj. 

Our final RFI mitigation technique is manual - in 
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our results it is easy to see high-cr hits that are a result 
of RFI: these hits appear as simultaneous detections at 
many dispersion measures. 



3.3. Detection of Giant Pulses from the Crab Nebula 

A suitable test of transient detection capability is 
to observe the Crab pulsar and attempt to detect giant 
pulses from it. 

Figure |5] shows a diagnostic plot generated from a 
one-hour Crab observation. The data is an incoherently 
summed set from the 35 best inputs. The diagnostic 
plot was generated after the raw data had been dedis- 
persed using a range of dispersion measures from 5 
to 200. The Crab pulsar has dispersion measure » 57 
cm~^ pc, so three giant pulses from the Crab pulsar 
can be easily identified in the lower plot. The giant 
pulses appear only at the expected dispersion measure, 
whereas wideband RFI appears across a wide range of 
dispersion measures. 
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Figure 5: Diagnostics on data taken from the Crab pulsar in a 60- 
minute observation conducted on 22 December 2007. Tlie data was 
dedispersed using dispersion measures ranging from to 200 cm"^ 
pc. Top-left: single-pulse SNR liistogram. Top-centre: noise appears 
at all DMs, but bright pulses (SNR> 7) from the Crab correctly appear 
at DM X 57 cm"'' pc. Top-right: inset of Figure|6] Bottom: pulse de- 
tections plotted on the DM versus time plane. Higher SNR detections 
appear as larger circles. Three giant pulses from the Crab are clearly 
visible in this plot. 

A frequency vs. time plot of the raw (summed) data 
at the time when the brightest giant pulses was detected 
is shown in Figure |6] The pulse is clearly visible in 
the data, and a fit to the dispersion measure shows that 
the pulse is, nearly without doubt, from the Crab (as 
opposed to RFI). 



Crab Pulsar: Giant Pulse Deleclion (sigma = 15.75) 
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Figure 6: A giant pulse from the pulsar in the Crab nebula. This pulse 
was detected in a 60-minute dataset taken on 22 December 2007. The 
dispersion of the pulse, correctly con'esponding to DM = 56.78 cm"' 
pc, is clearly visible. 

4. Astropulse 

4.1. Algorithm 

Astropulse uses a unique approach to the problem of 
searching for dispersed pulses. This problem is well 
known to pulsar astronomers, and many solutions have 
been devised over the years. The traditional method 
of searching for pulses of unknown dispersion measure 
(DM) is through incoherent de-dispersion, which is ba- 
sically a filter bank - for each channel, the power is mea- 
sured and integrated, then appropriate delays are added 
to compensate for the dispersion. This can be made very 
efficient through the use of a binary tree type algorithm. 
We use a coherent de-dispersion algorithm, which pre- 
serves the phase of the signal, thus increasing sensitiv- 
ity. It also allows us to have much better time resolu- 
tion (down to the band limit at 0.4 fis ) than incoherent 
methods, whose theoretical maximum time resolution 
is determined by the individual channel bandwidths. (In 
practice, the time resolution for pulsar searches is typi- 
cally 100 fis or greater.) The problem with coherent de- 
dispersion is that it is very computationally intensive. 
We are able to afford this by implementing the pulse 
search in a distributed computing environment. 

The coherent de-dispersion process is basically a con- 
volution. The time domain data needs to be convolved 
with an appropriate chirp function in order to remove 
the dispersion smearing effect. 

t 

fch{t) = exp(2;riv(f)f) = exp(2;ri — Avf ) 

of 

Here, 5t is the amount of time stretching caused by 
dispersion over a bandwidth of Av, which can be deter- 
mined from equation [T| 
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The most efficient way to do this is through the use 
of FFT convolution. Here is a brief description of the 
Astropulse detection algorithm: 

1. Take 13.4s time-domain data, perform FFT on 
32k-sample chunks. Overlap chunks by 50% to re- 
cover pulses which span a chunk boundary. Save 
results on disk for later use. 

2. For each DM, 

(a) Multiply frequency domain data by correct 
chirp. 

(b) Inverse FFT to go back to the time domain. 

(c) Threshold the (de-dispersed) time data, 
recording strong pulses. 

(d) Co-add adjacent bins, and threshold to look 
for broad pulses. 

(e) Search repeating pulses by using a folding or 
harmonic search algorithm. 

3. Repeat for next DM. 

4. When DM range is covered, start at step 1 for next 
13.4 s of data. 

The optimal step size for DM is to step by the 
amount of dispersion which would cause 1 extra time 
sample length of stretching. Using the above ex- 
pression for 6t - 0.4 yKS , this gives us a DM res- 
olution of 0.055 pc cm"^ . We want to cover DMs 
from 55 pc cm"^ to 830 pc cm"^ , which corresponds to 
14,000 DM steps at this resolution. (DMs smaller than 
about 55 are inaccessible to us, or to any search that 
employs one-bit sampling.) 

The maximum DM we wish to search determines 
the length of the FFTs. We want to look at DM up 
to 830 pc cm"^ , which will cover most of the galaxy. 
Using the same formula, this is a time stretch of St - 
5600 jUS , or about 14,000 samples. Our FFT size should 
be at least twice this big, so we will use 32k point FFTs. 

A FFT takes about 5A^ log floating point operations 
(FLOPs) to compute, where is the number of points 
in the FFT. We are doing = 2'^ point FFTs, and we 
need to do one inverse transform for each DM, plus one 
original forward transform. In addition, there are an- 
other multiplies per DM for the chirp function, and 
a factor of 2 for the 50% overlap. This makes our total 
number of operations for a 32k sample chunk of data: 

Nops = 2{Ndm + \)N\ogN + 2NdmN = 1.8 x 10'° 

To calculate how much computation we would need 
to analyze the data in real time, divide the above num- 
ber by 0.013 s (length of 32k samples). We need an- 
other factor of 14 from multiple beams (7) and multiple 



polarizations (2). Inefficiencies in memory allocation 
processes may contribute another factor of 5. In all, 100 
TeraFLOP/s or more may be required. For comparison, 
an average desktop PC can compute at a speed of about 
1 GigaFLOP/s. 

This analysis does not take into account the amount 
of computation required for folding. Astropulse folds 
on two different time scales, spending enough time on 
these to triple the overall computation requirement. 

This computation would take far too long on one 
computer, or even on a modest sized Beowulf cluster. 
Our solution to the problem is to distribute the data 
publicly to volunteers who download a program which 
will analyze it on their home computers. This is set 
up as a screen saver, which will turn on and start com- 
puting when their computers would otherwise be idle. 
(|6l This approach has been remarkably successful for 
SETI@home, which has an average computation rate of 
about 40 TeraFLOP/s. 

We are able to implement this though the BOINC 
software platform, also developed by our group at 
Berkeley. The package, the Berkeley Open Infras- 
tructure for Network Computing (BOINC), is a gen- 
eral purpose distributed computing framework. It takes 
care of all the "management" aspects inherent to a dis- 
tributed computing project - keeping records of user ac- 
counts, distribution of data and collection of results over 
the network, error checking via redundant processing, 
sending out updates of the science code, etc. Many 
lessons learned over the course of running the origi- 
nal SETI@home project went into BOINC, so that new 
groups wishing to start computing projects don't have 
to "re-learn" these. 

An interesting aspect of running a public distributed 
computing project is that it has education and pub- 
lic outreach automatically built into it. People invest 
their computer time into the project, and as a result 
become interested in learning more about the science. 
SETI@home has been enormously successful in this re- 
gard - the website averages 1.25 million hits per day, 
and the project has spawned many internet discussion 
groups and independent websites. It is also used as part 
of the science curriculum by thousands of K-12 teachers 
nationwide. 
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